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Abstract 

In this work we examine the operation of a node relaying packets from a number of users to a destination 
node. We assume multi-packet reception capabilities for the relay and the destination node. The relay node can 
transmit and receive at the same time, so the problem of self interference arises. The relay does not have packets 
of its own and the traffic at the source nodes is considered saturated. The relay node stores a source packet that 
it receives successfully in its queue when the transmission to the destination node has failed. We obtain analytical 

expressions for the characteristics of the relay's queue (such as arrival and service rate of the relay's queue), the 

stability condition and the average length of the queue as functions of the probabilities of transmissions, the self 
interference coefficient and the outage probabilities of the links. We study the impact of the relay node and the self 
^ interference coefficient on the throughput per user-source as well as the aggr-egate throughput. 

m 

>- I. Introduction 

(N 

The classical relay channel was originally introduced by van der Meulen yj . The first works on the relay channel 
^ were based on information theoretical formulations as in ||2| and |[3|. Recently several works investigated relaying 
capabiUty at the MAC layer ||3|, Q, ||5|, |[6|. The classical analysis of random multiple access schemes like slotted 
O ALOHA ||7| has focused on the so called collision model. Random access with multi-packet reception (MPR) has 
^ attracted attention recently |[8|, Q, 1 10 1. All these previous approaches come together in the model that we consider. 



^ In wireless networks when a node transmits and receives simultaneously the problem of self interference arises. 
^ Information theoretic aspects of this problem can be found at the work of Shannon on pi[ , although the capacity 
^ region of the two-way channel is not known for the general case 1 12 1. There are some techniques that allow the pos- 



sibility of perfect self interference cancelation 1 12 1. In practice though, there are technological limitations 1 13 1- 1 14 1 



which can limit the accuracy of the self interference cancelation. Various methods for performing self interference 



cancelation at the nodes' receivers can be found in [ 15] and | [T6| . The conclusion is that there is a trade off between 
transceiver complexity and the accuracy of the self interference cancelation. 

In this work we examine the operation of a node relaying packets from a number of users-sources to a destination 



node as shown in Fig. [TJ and is an extension of |17| and |18|. We assume MPR capability for the relay and the 
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destination node. The relay node can transmit and receive at the same time. We assume random access to the 
channel, time is considered slotted, and each packet transmission takes one time slot. The wireless channel between 
the nodes in the network is modeled by a Rayleigh narrowband flat-fading channel with additive Gaussian noise. A 
user's transmission is successful if the received signal to interference plus noise ratio (SINK) is above a threshold 
7. We also assume that acknowledgements (ACKs) are instantaneous and error free. The relay does not have packets 
of its own and the sources are considered saturated with unlimited amount of traffic. We do not consider any specific 
self interference cancelation mechanism, because it is out of the scope of this work. The self interference cancelation 
at the relay is modeled as a variable power gain. 

We obtain analytical expressions for the characteristics of the relay's queue (such as arrival and service rates), we 
study the stability condition and the average length of the queue as functions of the probabilities of transmission, 
the self interference coefficient and the outage probabilities of the links. We study the impact of the relay node and 
the self interference coefficient on the throughput per user-source and the aggregate throughput. 



Section [11] describes the system model, in Section [III] we study the characteristics of the relay's queue and we 
derive the equations for the throughput per user and the aggregate throughout. We present the numerical results in 
Section IV and, finally, our conclusions are given in Section |V] 



II. System Model 

A. Network Model 

We consider a network with sources, one relay node and a single destination node. The sources transmit 
packets to the destination with the cooperation of the relay; the case of A^ = 2 is depicted in Fig. [T] We assume 
that the queues of the two sources are saturated (i.e. there are no external arrivals but unlimited packet volume in 
the buffers); the relay does not have packets of its own, and just forwards the packets that it has received from 
the two users. The relay node stores a source packet that it receives successfully in its queue when the direct 
transmission to the destination node has failed. We assume random access to the channel. Each of the receivers 
(relay and destination) is equipped with multiuser detectors, so that they may decode packets successfully from 
more than one transmitter at a time. The relay node can receive and transmit packets simultaneously. 



l-p. 





Fig. 1. The simple network model 



Fig. 2. Markov Chain model for the two- 
user case 



B. Physical Layer Model 

The MPR channel model used in this paper is a generalized form of the packet erasure model. We assume that 
a packet transmitted by i is successfully received by j if and only if SIN R{i, j) > 7^, where 7^ is a threshold 
characteristic of node j. The wireless channel is subject to fading; let Ptx{i) be the transmitting power at node i 
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and r{i,j) be the distance between i and j. The received power by j when i transmits is Prxihj) = j)g{i, j) 
where A{i,j) is a random variable representing channel fading. Under Rayleigh fading, it is known | [l9[ that A{i,j) 
is exponentially distributed. The received power factor g{i,j) is given by g{i,j) = Ptxii){f{i, j))'"^ where a is 
the path loss exponent with typical values between 2 and 4. We model the self interference by a scalar g G [0, 1]. 
We refer to the g as the self interference coefficient. When g = I, no self interference cancelation technique is used 
and g = when there is perfect self interference cancelation. The success probability in the link ij is given by: 

f ^j^j ^n . f -wc \-rn rr f, , v{k,j)g{k,j)\''^ 
''^ V v{i,])g{i,])J i.cT\7- V v{i,j)g{i,j) J 

where T is the set of transmitting nodes at the same time, v{i,j) is the parameter of the Rayleigh random variable 
for fading; m = 1 when j ^ T and m = else. The analytical derivation for this success probability can be found 
in 120) and ||2TJ. 

III. Analysis 

In this section we derive the equations for the characteristics of the relay's queue, such as the arrival and service 
rates, the stability conditions, and the average queue length. We will provide an analysis for two cases: first, when 
the network consists of two users (non-symmetric) and the second is for n > 2 symmetric users. 



A. Two-user case 

I) Computation of the average arrival and service rate: The service rate is given by: 

= go(l - - q2)PQ/o + qoqi{l - q2)Po/o,i + QoQ2{l - qi)Po/o,2 + 909192^0/0,1,2 (2) 

where qq is the transmission probability of the relay given that it has packets in its queue, qi for i 7^ is the 
transmission probability for the i-th user. The term P^y . ^ is the success probability of Unk ij when the transmitting 
nodes are i and k and can be calculated based on ([T]). 
The average arrival rate A of the queue is given by: 

A = P(Q = 0)Ao + P(Q>0)Ai (3) 

Where Aq is the average arrival rate at the relay's queue when the queue is empty and Ai when it's not. 
-^0 = ^1 + 2r2, where is the probability of receiving i packets given that the queue is empty. The expressions 
for are given by: 

r? = <?l(l - 92)(1 - Pi'/JP?/! + 92(1 - 5i)(l - P|/2)P°/2 + 9192(1 - P^/ 1, 2) P?/ 1^2^1,2+ 
+9192(1 - Pi/^,2)P2/l,2Pl/1.2 + 5192(1 - PlA,2)Pl/l,2(l - P2A,2)(1 " ^2/1,2) + '?1'?2(1 - P2A,2)^2/1,2 (1 " PlA,2)(l " P?/l,2) 

4 = (71(72(1 - Pi'/i,2)(l - ^2A,2)^'V1,2^2/1,2 (5) 

Accordingly, Ai = r} + 2r2, where rj is the probability of receiving i packets when the queue is not empty. The 
expressions for the rj are lengthy and given by: 
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rl = {l- qo)qi{l - q2)(l - Pt/i)P?/i + 9091(1 - <?2)(1 - P^/o,i)P?/OA + (1 " 90)92(1 - 9l){l - Pi/2)P2/2 + 9092(1 - 9l)(l - P2/o,2)P2/0,2 + 

+ (1 - 90)9192(1 - fvi,2)A°/1.2(l - ^2/l,2)(l - ^■2A,2) +909192(1 - A%,1,2)^'?/0,1,2(1 " ^2A),1,2)(1 ^ ^/0,1,2) + 

+ (1 - 90)9192(1 - P,%2)P?/1,2P2/1,2 + 101112(1 ~ /'i%,l,2)^r/0.1,2^2/0,l,2 + 

+ (1 - 90)9192(1 - P2A,2)^'2/1.2(1 - ■PlA,2)(l - ■PlA,2) + 909192(1 " ^'2%,1,2)^2%,1,2(1 " ^V0,l,2)(l - P°^o,l,2^ + 

+ (1 - 90)9192(1 - P^/l,2)P2/1.2Pl/l,2 + 909192(1 - ^'2%,1,2)^1/0,1,2-Pi%,l,2 

(6) 

4 = (1 - qo)qiq2{l - ^i"/i,2)^'iA,2(l " ^2/1,2)^2/1,2 + 909192(1 - ^'i%,i,2)^i%,i,2(l " ^2/0,l,2)^2%,l,2 (7) 

In Fig. |2] we present the discrete time Markov Chain (DTMC) that describes the queue evolution. Each state is 
denoted by an integer and represents the queue size at the relay node. The transition matrix of the above DTMC 
is is a lower Hessenberg matrix given by: 

^ ao bo 
ai bi bo 
02 62 ^1 
63 62 
63 

Where oq = 1 — — P2, ^i = Pi, 02 = P2> ^0 = P-i and ftj+i = p} i = 0, 1, 2, 3. The quantity p^ (p^) is the 
probability that the queue size increases by i packets when the queue is empty (not empty). Note that p^ = r?, 
because when the queue is empty the probability of i packets arriving is the same with the probability that the queue 
size increases by i packets; when the queue is not empty however, this is not true. For example the probability of 2 
packets arriving is not the same with the probability of increasing the queue size by 2; this is because both arrivals 
and departures can occur at the same time. The expressions for the pj are also given by lengthy expressions listed 
below: 

pll = 90(1 - gi)(l - 92)^0% + 90(1 - 91)92 P(f/0,2-P2/0,2 + 90(1 - 9l)92P(f/o,2(l " ■P2/0,2)(l " ■P2/0,2) + 
909l(l - 92)Po'/0,l^lVl + lona - 92)Po''/0,l(l - Pi/OaK-^ - P?/0,l) + 9o9ig2PoVl,2-P'/0,1.2^2'/0.1.2+ 
+909l92Po'/0,l,2(l - ■PV0,1.2)(1 - Pi%,i,2)(l - ^2/0.1,2)(l " ■P2%,1,2) + 909l92P(5^o,l,2^'vo,l,2(l " ■P'/0,l,2)(l " ^2%,1,2) + 

909192 Pq'/o,!, 2(1 - ■PV0,1,2)(1 - Pi%,i,2)-P2/0,1,2 

pI = i-pIi-p\-pI (10) 

pi = (1 - 90)91(1 - 92)(1 - Pi'*/i)P{'/i + (1 - 90)9192(1 - Pvi,2)Pl/1.2^2A.2 + 
(1 - 90)9192(1 - P{'/i,2)Pi°/l,2(l - ^2A,2)(1 - ^2/1,2) + (1 - 90)(1 - 91)92(1 - P2a)^2/2 + 
(1 - 90)9192(1 - P2A,2)^2/1.2Aa,2 + (1 " 9o)9l92(l - P2A,2)^'2/1,2(1 ' ^'lA,2)(l " ^"71,2) + 
+909192Po%,i,2(l - n'/0,l,2)^V0.1,2(l - ^2%,1,2)^2%,1,2 + 9091(1 " 92)(1 " Po%,i)(l " P,%i)P?/o,l+ dD 
+909192(1 - Po%,l,2)(l - ^'V0.1,2)^V0.1,2^2/0.1,2 + 909192(1 - Po%,l.2)(l " n%,1.2)^l/0,l,2(l " ■P2/0,1.2)(l " ^2%,1.2) + 
+9092(1 - 9l)(l - Po%,2)(l - ^2A),2)^2%,2 + 909l92(l - PoVl,2)(l " ^2'/0,l,2)-P2%,l,2^V0,1.2 + 
+909192(1 - Po%.l,2)(l - ^^2%,1.2)^2%,1,2(1 " Pt/0,l,2)i^ " ■Pi°/0.1,2) 

P2 = (1 - 90)9192(1 - PlA,2)Pl/1.2(l - ^2A.2)-P2A.2 + ^09192(1 " Po/o,l.2)(l " ^'vo,l,2)-Pr/0,l,2(l " ■P2%,l,2)-P2/0.1,2 ('2) 






^0 

bi 
b2 



(8) 
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The difference equations that govern the evolution of the states are given by: 

Ps = s ^ Si = aiSQ + ^ bi^j+iSj (13) 
i=i 

We apply the Z-transform technique to compute the steady state distribution, i.e. we let 

2 3 oo 

A{z) = ^ a^z-\B{z) = ^ biZ-\ S{z) = ^ s^z"* (14) 

i=0 i=0 i=0 



It is known that |22|: 



z-^Aiz)-Biz) 

Siz) = so (15) 



It is also known that the probability of the queue in the relay is empty is given by p2 | 

1 + B'{1)-A'{1) 



^'(Q = 0) = :^4S#H77T (16) 



The expressions of A' (l) and -B (1) are: 

(2 \ ' 2 2 2 

E = " ^ ia.z-(^+i) ^A'{l) = -Y,m^A'{l) = -Y,w' = -^o (17) 
i=0 / i=l i=l i=l 

/ 3 \ ' 3 3 

^'(^) = E = " E ^^^^"^'^'^ ^ ^' (1) = - E = - 262 - 363 = -1 + P-i -p\- 2pl 

\i=0 / 2=0 i=0 

(18) 

Then the the probability of the queue in the relay is empty is 

PiQ = 0)= (19) 

pil - Pi- 2p^ + Ao 

So, the average arrival rate A is given by: 

A= y\ Ao + - /° , , Ai (20) 

2) Condition for the stability of the queue: An important tool to determine stability is Loyne's criterion (231, 
which states that if the arrival and service processes of a queue are jointly strictly stationary and ergodic, the queue 
is stable if and only if the average arrival rate is strictly less than the average service rate. If the queue is stable, 
the departure rate (throughput) is equal to the arrival rate. \i < r\ + 2r\ < /x where r\ = {1 — qo)Ai + qoBi, 

= {1 — qo)A2 + qoB2 and /x = q^A. The expressions for A, Ai, Bi are given by: 

= - g2)(l - Pt/^)Pl/^ + '72(1 - gi)(l - Pii^)P°/^ + 9192(1 - Pt,^,2)Pin,2^^ " ■P2/i,2)(l " ^2A,2) + 

+gi'72(l - Pi/^,2)PinA,2 + 9192(1 - ^'2A.2)^2A,2(1 " Aa,2)(1 " ^1/1,2) + 9l92(l " P2A,2)^2/l,2^Vl.2 

(21) 

Bl = gi(l - 92)(1 - + <72(1 - gi)(l - P2%,2)^2%.2 + 9192(1 " P^ / ,1 ,2^ ^1 / ,1 ,2^^ " ■P2'/0,l,2)(l " ^'2%,1,2) + 

+9192(1 - nVl,2)n%,l,2^2%,l,2 + 9192(1 - ^'2'/0,l,2)^2%,l,2(l " ■Pv0,l,2)(l - P{'/o,i,2) + 

+9192(1 - ■P2/0,l,2)-^2/0, 1,2-^1/0, 1,2 
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B2 = gig2(l - -Pf/o, 1,2) ^1/0,1,2(1 - ^2%,1,2)^2%,1,2 



^ = (1 - gi)(l - q2)P^/o + - Q2)Po/0,l + 92(1 - gi)^o/0,2 + '?l'?2</o,l,2 

Then the values of qo for which the queue is stable is given by qomin < qo < ^, where: 

Ai + 2A2 



qon 



A + A1+2A2-B1- 2B2 



3) Average queue size: The average queue size is given by |22|: Q = —5 (1) where 5 (1) 
expressions for K{z) and L{z) are given by: 



(22) 
(23) 

(24) 



sojjr^. The 



K{z) = [-z-^A{z) + z-^A{z) - B\z)^ (z-i - B{z)) - {z-^A{z) - B{z)) (-^-^ _ b\z)^ (25) 



L{z) = {z-' - B{z)f 

Then i<r"(l) and are given by: 

^ = f2yl(l) - 2A'(1) + A"(l) - i?"(l)) f-1 - B\l)\ -[2-B''{l)]{ -A{1) + A [l] - B' {I) 



(26) 



L"(z) = [2 (z-i - B(z)) _ ^ L"(1) = 2 (-1 - B\l) 

The values of A (1) and B (1) are: 



(27) 
(28) 



E 

j=i 



laiZ 



^ i(i + l)aiZ-(*+2) ^ A" (1) = 2p? + 



The average queue size is given by: 

Q 



(29) 



(z) = ( - E =Y.^{i + l)biz-^'^^^ ^B"{1) = 2- 2p\^ + Ap\ + lOp^ (30) 

V i=l J 1=1 



(31) 



2(pi + 2pl - pl^)ipl, -p\- 2p\ + Ao) 

4) The throughput per user and the aggregate throughput: The throughput rates /xi,/i2 for the users 1,2 are 
given by: 



IJ^l = qoP {Q > 0) qi{l - q2) + (1 - Pt/OA)Pl/oA + (Q > \Pt/ox2 + " n%,l,2)n%,l,: 



+ 



+ [1 - qoP {Q > 0)] qiil - q2) [Pf/i + (1 - Pt/i)Pi/i\ + [1 " QoP {Q > 0)] qiq2 [Pt/1,2 + (1 " P'l/l,2)P'l/l,2 

(32) 
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M2 = g0^'(g>0)(?2(l-(?l) P2/0,2 + i^-P2/0,2)P2/0,2 + QoP (Q > 0) 



pd I (-] _ pd \ pO 

-^2/0,1,2 ^ V-L -^2/0,1,2^^2/0,1,2 



+ 



+ [1 - qoP {Q > 0)] (?2(1 - qi) [P^/2 + (1 - -P2'/2)^2°/2] + [1 - QOP {Q > 0)] 9192 [Pi/i,2 + (1 " -P2/l,2)^'2/l,2 

(33) 

In the equations above we assume that the queue is stable, hence the arrival rate from each user to the queue is a 
contribution to its overall throughput. The aggregate throughput is utotai = )Ui + A*2. Notice that the throughput per 
user is independent of qo as long as it is in the stability region. This is explained because the product qoP {Q > 0) 
is constant. The proof is straightforward and thus is omitted. 

When the queue is unstable however, the aggregate throughput is the summation of all the direct throughput 
between the users and the destination plus the service rate of the relay. 

B. N-symmetric users 

We now generalize the above for the case of a symmetric n-users network. Each user attempts to transmit 
in a slot with probability q; the success probability to the relay and the destination when i nodes transmit are 
given by Po,i, Pd,i respectively. There are two cases for the Pd,i , Pd,i,o^ Pd,i,i denoting success probability when 
relay remains silent or transmits respectively. The above success probabilities for the symmetric case are given 

by Pd,,j = Pd {r^y {r^y , j = 0, l and /3 = ^ > l. Pod, = Pod [^^J , Po = cxp (-^), 

Pd = exp (^-^), Pod = exp There are two cases for the Po,i, -Po,i,o, -Po,j,i denoting success probability 

when relay remains silent or transmits respectively. The success probabilities are given by Po,i,o = Po (^j^^ 

and Po,i,i = (1 + 7o^o5)~^ (tTto) where vq is the distance between the users and the relay, a is the path 
loss exponent and g is the self interference coefficient. 

1 ) Computation of the average arrival and service rate: The service rate is given by the following equation: 

^ = E [f) ^o^'(l - 'ir''Pod,k (34) 
The average arrival rate A of the queue is given by: 

A = P (g = 0) Ao + P (g > 0) Ai (35) 

Ao = J2k=i ^''"k where the is the probability that the relay received k packets when the queue is empty, the 
expression for is given by: 

^k = J2 (") (l) 9^(1 - 9r"Po,,o (1 - Pd,i,o)'' [1 - Po,.,o(l - Pd,i,o)Y-\ l<k<n (36) 

Ai = Ylk=i ^^fe where the is the probability that the relay received k packets when the queue is not empty 
and is given by: 



n 



i=k 



tI = {1- go)E (^) [Mil - qr-'P^,i,o (1 - Pd,i,o)' [1 - PoM^ - Pd,i,o)r'+ 



+50 E Q Qq\l - Qr-%'i,i (1 - Pd,i,i)' [1 - Po,iA^ - Pd,i,i)r^ l<k<n 



(37) 



The elements of the transition matrix are given by: = p^, bo = pL^, bi = Pq and bk+i = p\yk > where: 



i=k 



n\ I 



i j \ k 



q\l - qr-'P^,,o (1 - Pd^flT [1 - ^'o,,o(l - Pd,fi)r\ l<k<n 



i—k 



pI = (1 - ^o) E h) (l) q\l - qr-'Po%,o (1 - ^^^,.0)' [1 - Po.,o(l - Pd,,o)Y-''+ 



(38) 



(39) 



-90 



E 

i=k 



i=k+l 



i J \k 



q\l - qr-\l - Pod,)PLl (1 - Pd,,l) [1 - ^'o,,l(l - Pd^)] + 



+^0 E [i)[k + 1 j'^'^^ " «)""^0'^-^^o,S (1 - ^'^.u)'^' [1 - ^0,i,l(l - PdM)] 



i-k-l 



(40) 



i=l 



(41) 



The probability that the queue in the relay is empty is given by ( 16 1, the expressions for A (1) and B (1) are: 



A(z) 



^1=0 



ajz 



E 



-(i+l) 



A'il) = -^ia,^ A'{1) = -^ip^ = -Ao (42) 



i=l 



/n+1 



i?'(z)= 



n+l 



n+1 



n+l 



^ 5' (1) = - ^ i6, = - ^ = -1 + _ ^ i;pi (43) 



\i=0 / i=i i=i 1=2 

Then the probability that the queue in the relay is empty is given by: 



1=1 



P-l - E ^Pi 



P(Q = 0) 



i=le 



P-1 -JZ'^Pi + 

j=l 



(44) 



2) Condition for the stability of the queue: Xi < ^ ^ X]fc=i ^''fc ^ ^ where = (1 — qo)Ak + (/o-B^ and 
^ = q^A. The expressions for A, A^, are : 



i=k 



i ) \k 



q\l - qr-'Po% (1 - ^d,.o)'^ [1 - ^'o,,o(l - Pd,fl)] 



i—k 



(45) 



i=k 



A = Y.^I^^^'i^-ir-'Pod,k 



(46) 



(47) 
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The values of go for which the queue is stable is given by gomm < 9o < 1> where: 



qor, 



k=l 



A + "^kAk-Yj 



(48) 



k=l 



k=l 



3) Average queue size: As we showed in the Section 



III-A3 



the average queue size is given by: Q = —S (1) 



where 5 (1) = sq ^„ i^) . The expressions for K" (1) and L" (1) are given by (27i and (28 1. The expressions for 



A"{1) and ^"(1) are 
A"(z) = 



+ l)aiZ-('+') ^ X' (1) = Y,i{i + l)ai = 5^ i(i + 



(49) 



i=l 



1=1 



i=l 



i=l 



B (z) 



n+1 \ n+1 



n.+l 



J2 + l)biZ-'^'+^^ B" (1) = ^ i(i + l)hi = 2 - 2pii + Y,^{i + ^)p\ 

(50) 



1=1 



i=l 



i=l 



Following the same methodology as in Section III-A3 we obtain that the average queue size is given by: 

/ n \ ( ^ \ 



1=1 



i=l 



(51) 

2fE^^^-^-il ki-E^^^' + ^oj 

4) The throughput per user and the aggregate throughput: The throughput per user for the network with the 
relay when the queue is stable is given by: 



n-l / _ i\ 
I — n \ / 



\n—l—k 



Pr 



d,fc+l,l 



(1 — -P(i,A:+l,l)-fo,A;+l,l] + 



fc=0 
n-l 



(52) 



+ [1 - goP (Q > 0)] E ( / ) - 9)"-'-' [Pd,k+ifl + (1 - Pd,fc+i,o) Po,fc+i,o] 

fc=0 ^ 



When the queue is unstable though, the throughput per user is given by the summation of the direct to the 
destination throughput plus the service rate of the relay divided by the number of the users n. The aggregate 
throughput is p.totai = nfi. 

IV. Numerical Results 

In this section we present numerical results for the analysis presented above. To simplify the presentation we 
consider the case where all the users have the same link characteristics and transmission probabilities. The parameters 
used in the numerical results are as follows. The distances in meters are given by = 130, tq = 60 and tq^ = 80. 
The path loss is a = 4 and the receiver noise power rj = 10^^^. The transmit power for the relay is PtxiO) = W 
mW and for the i-th user Ptxii) = 1 mW. We used 7 < 1 because it is possible for two or more users to transmit 
successfully at the same time. 



The figures 3(a) and 5(a) present the throughput per user versus g (the self-interference coefficient) for various 



values of 5,7 and n. The figures 3(b) and 5(b) show the aggregate throughput versus g 
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(a) Throughput per user vs self interference coefficient (b) Aggregate throughput vs self interference coeffi- 
cient 

Fig. 3. Throughput per user and aggregate throughput vs the self interference coefficient for 7 = 0.2 
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(a) Throughput per user vs the number of the users 
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(b) Aggregate throughput vs the number of the users 



Fig. 4. Throughput per user and aggregate throughput vs the number of the users for 7 = 0.2 





(a) Throughput per user vs self interference coefficient (b) Aggregate throughput vs self interference coeffi- 
cient 

Fig. 5. Throughput per user and aggregate throughput vs the self interference coefficient for 7 = 0.6 



The figures 4(a) and 6(a) present the throughput per user versus the number of the users in the network for 



various values of g,and for 7 = 0.2 and 7 = 0.6 respectively. The figures 4(b) and 6(b) show the aggregate 
throughput versus the number of the users. When 7 = 0.2 we observe that for g = 10^^° and g = 10^® (almost 
perfect self-interference cancelation) the relay's queue is unstable for relative small number of users. The previous 
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result is because the small value of 7 is more likely more transmissions from the users to the relay to be successful, 
but at the same time the relay can transmit at most one packet per time slot. For 7 = 0.6 the queue is never unstable 
for the parameters described in the figures, and for g = 10^^'^ and g = 10~^ the advantages in term of throughput 
are obvious compared to no self interference cancelation. 

The Fig. |7(a) and Fig. 7(b) show the gomm vs n for 7 = 0.2 and 7 = 0.6 respectively. Note that gomm < g < 1, 
so when qomin > 1 then the queue is unstable as in case of 7 = 0.2 for g = 10"^° and g = 10"^. 




(a) Stability threshold vs the number of the users for (b) Stability threshold vs the number of the users for 
7 = 0.2 7 = 0.6 

Fig. 7. Stability threshold vs the number of the users 



V. Conclusions 

In this paper, we examined the operation of a node relaying packets from a number of users to a common 
destination node. We assumed MPR capability for the relay and for the destination node. We studied a multiple 
capture model, where a user's transmission is successful if the received SINR is above a threshold 7. The relay 
node can also receive and transmit simultaneously, so the problem of self interference arises. 

We obtained analytical expressions for the relay's queue characteristics such as the stability condition, the values of 
the arrival and service rates, the average queue size. We studied the throughput per user and the aggregate throughput, 
and found that, under stability conditions, the throughput per user does not depend on the relay probability of 
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transmission. We studied the impact of self interference coefficient on the throughput per user and the aggregate 
throughput of the network. 

We showed that for perfect self-interference cancelation, the advantages are obvious. Another interesting result 
is that the self interference coefficient plays a crucial role when 7 is small (and g tends to zero) because it can 
easily cause an unstable queue. 

Future extensions of this work should include users with non-saturated queues i.e. sources with external random 
arrivals, a relay node with its own packets and different priorities for the users. 
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